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TECHNICAL MEMORANDUM X-64851 


DESIGN OF A DIGITAL CONTROLLER FOR SPINNING FLEXIBLE SPACECRAFT 

Summary 

A new approach to digital control system design is applied to the analysis and de- 
sign of a practical onboard digital attitude control system for a class of spinning vehicles 
characterized by a rigid body and two connected flexible appendages. The approach used 
is to design a continuous -data control system that will provide a satisfactory system re- 
sponse. Then, using the digital redesign method, a digital controller with onboard digital 
computer is designed to provide a digital control system whose states are simitar to those 
of the continuous system at sampling instants. The simplicity of application of this ap- 
proach is indicated by example. The example, using spinning Skylab parameters, is 
used to substantiate the conclusions. 

I. Introduction 

In 19 70, NASA initiated a set of studies of several Skylab configurations, includ- 
ing a spinning vehicle configuration. (The spin would provide an artificial-gravity 
environment to assess and compare the physiological benefits and problems of prolonged 
zero-gravity and artificial-gravity environments.) The study proceeded generally along two 
mutually reinforcing directions. In one case, a large-scale digital simulation of the 

Presented as Paper 73-894 at the AIAA Guidance and Control Conference, 

Key Biscayne, Fla., Aug. 20-22, 1973. 

Indexes Categories: Spacecraft Attitude Dynamics and Control; Navigation, 

Control, and Guidance Theory. 
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spinning system of connected flexible and rigid structures was developed, first devis- 
ing a program for obtaining the eigenvalues and eigenvectors and then simulating the 
truncated equations for motion* (1) 9 (2) 9 (3) Ib order to develop faith in the response 
of the large scale simulation and to obtain characteristic numerical values for the 
control system parameters , a simplified analytically-tractable model was developed and 
analyzed in detail. (4), (5) While the model was developed to meet Skylab demands, it is 
applicable to other spinning bodies with flexible attached members. It consists of a 
single rigid-core body with two attached flexible massless booms having tip masses 
(Fig* 1-1) * Active closed loop attitude control is implemented by applying control 
torques to the vehicle to compensate for the effect of disturbance torques. (In the Skylab 
example, existing Skylab onboard sensors and momentum exchange systems are used for 
generating control error signals and torques. ){6) 

The approach used is to design a continuous-data control system that wUl 
provide a satisfactory system response. Then using several recently developed 
techniques, the continuous control system n odel is redesigned to provide a digital 
control system. The objective of the digital redesign is to provide a system whose 
states are identical to those of the continuous system at sampling instants. 

H. Equations of Motion 

The spinning flexible vehicle is assumed to be modeled as shown in Fig. 1-1. 
For the entire vehicle, equations of displacement and rotation are written using Newton- 
Euler relations. The equations are based on the assumption that the axis of maximum 
moment of inertia of the overall vehicle does not coincide with the desired spin axis — 
usually an axis normal to the solar panels and pointed toward the sun. In that case, the 
overall mass distribution of the system must be altered to align the two axes-hence, the 
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inclusion of booms (assumed to be massless and flexible) with tip masses, ft is assumed 
that the booms are attached to the core rigid body at the axis of rotation of the overall 
system. The basic coordinate system's origin is at the steady-state mass center of the 
overall vehicle, and the axes coincide with the principal axes of inertia of the steady-state 

configuration of the vehicle. (In steady-state, the vehicle rotates about the major 
principal axis.) 

The nonlinear equation defining the rotational motion of the overall vehicle 
(Fig. l-l) is obtained from the Euler relation and becomes 

• ^ Q d 

T“ H = E . <^ + i> £ + yj x lt« u) + — J*_£ x £ dm (2-1) 

ere T is the applied torque vector, H is the angular momentum vector, i is the 
inertia diadic of the entire vehicle, w is the rotation vector of the entire vehicle, p is 
a generic vector joining the spinning steady-state location of the vehicle center of mass (CM) 
to each of the mass elements dm of the flexible appendages, the closed overdot ( • ) 
represents the time derivative in the inertial frame , and the open overdot ( «) represents 
the time derivative in a frame that is fixed with respect to the core rigid body. The 
nonlinear equation representing translational motion of each sth mass element (m S ) of 
the flexible appendages is obtained from the Newton relation and becomes 
F S = m S (X+c +_r S + u S ) 

-m S {X + c°+2ujxc+Sxc+cux(ajxc) + £ S + u° 8 

+ 2u)x(r S + uVj x (r S + ^ + u>x [<jjx (r®+ u 3 )]} (2-2) 

where F is the applied and structural connection force vector, X is the inertial 
position vector from an inertially fixed reference point to the instantaneous location of 
the CM, c is the vector from the CM to the nominal spinning steady state location of 
the CM S , r is a vector from CM g to the spinning steady state location of m s , and u s 


3 



Afl A flAl VlVli V\rr n4- A n#Jr< nin T i-* _ A L 

vw^i.^4. liwiu w±^ 0 £Ji Lining, oicauj S Lcltti lUUUliUIl UI in 


to its actual instantaneous location. 


If now these equations are linearized as perturbations about the spinning steady 
state, perturbation equations (2-3) and (2-4) result. Only the rotational relations have 
been retained, and the effect of the center of mass motion with respect to the body is 
neglected. The contribution of mass center motion due to flexible appendage motion and 
orbital dynamics is ignored. The resulting three linearized equations of motion are 


written explicitely in matrix notation: 

T 

T=Iw + Q Iw + wl Q + Q I Q + mfJ (2r^u B E - ru B - u B r^) 0 + 

m^zr u E-u r - ru )Q+ mru + mQ ru (2-3) 

where I is the moment of inertia matrix of the entire vehicle, 0= [[\ , Q a , Q 3 ] T is the 
nominal spin rate, w -[Wj^ , w 3 , w 3 ] T is the variation of spin speed, E is an identity 
matrix, m is the mass of each flexible boom's tip mass, u B = u 1 - u 2 , where u 1 and u 3 
represent the displacement of the tip masses of each of the two flexible booms from their 
spinning steady state location to their instantaneous location, and r = [0,r 2 , r 3 ] T represents 
the distance of the tip mass center of mass fn n the vehicle center of mass. A set of 
three equations describes the skew symmetric mode of the attached booms and tip masses: 
imi + (D + 2mQ)u + (K B + mQQ +K B )u B -m'r w + m(n w + w Q)r = 0 (2-4) 

5 

where K® = diag [k x , k a , k 3 ] and K B represent the elastic and geometric stiffness, 

respectively, of the attached flexible booms (represented as massless cantilever beams 

with attached discrete tip masses), and D B = diag [d x , d £ , d 3 ] is the structural damping. 

Stretching of the booms due to centrifugal force is represented by K 6 . 

s 

Active closed loop control is implemented by applying control torques to the 
vehicle. It is assumed that the 3-axis of the body-fixed coordinate frame must be held 
fixed in inertial space. This is achieved by applying control torques to the vehicle to 
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compensate for the effect of disturbance torques. These control torques must in some 
manner depend on error signals that are proportional both to the angle between the 
3-axis and the inertial reference vector and to its time derivative. This information 
can be generated if the spacecraft hardware includes, for example, sun sensors and 
rate gyros and the sun is used in place of an inertial reference. The sun sensor measure 
the angular rotations ^ (as also indicated for the simplified model of Fig. 1-1), and 
rate gyros measure the angular velocities w x , w s . The control torques T x , T a may be 
provided by momentum exchange devices such as control moment gyroscopes and 
augmented if necessary by a mass expulsion system; for simplicity it generally is 
assumed that T a = 0. (It is observed that the existing Skylab, which is used as a model 
to obtain numerical results herein, has onboard sun sensors, rate gyros, CMG's, and 

a mass expulsion system. )(6) A more complete derivation of the equations of motion is 
included in Ref. 7. 

A Linear control postulate (termed Control Law 1) can be formulated as 



where, from kinematic relations, 

w i] r- 1 0 °j *0i | ' 0 Q 0! *0 1 ‘ 

w 8 - 0 -1 0 0 2 + -Q 0 0 ! 0 2 

w 3 J j_ 0 0 -lj 0 3 j [ 0 0 oj _ 03 _ 


(2-6) 


An estimate of numerical values for the control gains may be obtained by an analysis of 
a simplified model of the vehicle. In that case - as shown subsequently - continuous- 
data control system gains for the output feedback of Eq. (2-5) may be determined by 
application of D-decomposition and parameter mapping techniques. (8) Alternately, 
optimal gains (Control Law 2) may be determined for complete state feedback by a 
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linear quadratic loss program to minimize a performance index. (9) Then the continuous 
control system model gains are altered to provide the desired system dynamic response 
when it is modeled as a digital system. 

When Eqs.(2-3), (2-4), and (2-6) are grouped in the conventional manner, the 
following set of equations results, assuming the steady state spin Q is about the 3-axis. 
In this case K® = diag [mfip, 0, mft 3 ]. 

o 

Mv M + D^’ + K/£ = - Fv 

where ' refers to differentiation with respect to t = nt, 

T 

$2 ? 03 9 Ma 9 ^2 > M-3 ] ♦ 

£"tVx. v 2» v 3 ] T = [T 1 /I 1 Q 3 , TaAiO 8 , T 3 /I 3 tf ] T 

and matrices M, D, K, and F are given by Eq. (2-10) through (2-13). The symbols used 
are defined as 


(2-7) 

(2-8) 

(2-9) 


K i (h ~ h)Ai, K z - (I 3 - I,)/i 2 , y. = 2rnr 2 2 /I. 


= (ii - u 3 )/2 r 2 


A i -d i /mn,(T i 2 -k i /mn 2 , (i = i, 2, 3 ),* = r 3 /r 2 . 

In the steady state, the principal axes of inertia of the total vehicle coincide with the 
1, 2, 3 axes and the principal moments of inertia are I L# I 2 , I 3 , respectively, with 
r i <J 2 <r a- The stiffness of the nonrotating booms in the directions of the 1, 2, 3 
axes is represented by k. and the structural damping by d.. 


M = 


“1 

0 

-r i 
0 
0 

J 


0 

1+K, 

1-K 2 

0 

0 

0 


0 

0 

0 

1 

1 

0 


0 

-€?i 

0 

n 

i 


iri 

o 

o 

0 

0 


-r, 

o 

y i 
o 

0 


(2-10) 
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o -(1+Ki) 0 

1+Ki 0 0 

0 0 0 

ooo 

-2£ 0 0 

0 - 2 £ 2 


0 0 

0 2 |y, 0 

0 0 YiA 3 

-2y 3 0 

-2 0 

A 2 0 


0 


K- 


K, 

0 

-y i 
0 

0 

H 

1 
0 


K ‘<lrl> 0 


f 

0 

* 

0 

0 0 

1 0 


0 0 


0 

0 

0 

0 

0 

0 

1 


0 -fri -r i 

£y } o o 

o o y,(<r 3 2 +i) 

ooo 

Crf 0 0 


0 


crl - 1 0 


0 0 
0 0 
0 0 


IT 


If control is exerted only about the 1- and 3- axes, then 


and 


-M'V- 


1 0 0 0 1 o o] T 

0 0 0 0 0 0 

0 0 0 1 0 0 


- f x 0 0 

0 0 0 

_ 0 -(l-fttcKj 


0 

0 

u 


—i T 


5fi 

0 

0 


< 

0 

0 


( 2 - 11 ) 


(2-12) 


(2-13) 


(2-14) 


(2-15) 
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(2-16) 

f x = l/(l -y 3 - ^^c). 

(2-17) 

c - (1 “ K 2 )/(l + Kj) 

(2-18) 

Eq. (2-7) may be cast in state equation form: 


z/= Az + Bv 

(2-19) 

where 


z = [y ; y] T , 

(2-20) 

A-[ 0 ^ 1 

l-M _1 K ! -M-iD J 

>2-21) 

b = [o 36 ;-m- 1 f] 

(2-22) 


and E s is the 6x6 identity matrix and Ogg is the 3x6 null matrix. 


It has been shown that when r 3 is sufficiently small that it may be ignored, the 
equations of motion (2-7) uncouple into two sets of equations: one describing the wobble 
motion (described by p, 3 amd w^, w s or cpj, cp^) and one describing the in-plane motion 
(described by p. £ , and w 3 or cp 3 ), (3) The resulting equation for wobble dynamics can 
be written as 

x ' «= Ax + B u (2-23) 

where 

x “ [<A* ^ 21 ^ 3 * 0'i» <&> H-i ] T (2-24) 
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0 

0 

0 

l 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

1 

Zi+h 

0 


0 

1+K, 

■TiA, 

i 

i-r, 


i -n 

0 

-k 2 

0 

“(l-Kz) 

0 

0 

Kill 

0 

y,-(cr 2 +l) 

o 

1+K, 


yrx 


i-r. 


1-Vi 

i-n J 


For control torque exerted only about the 1-axis, 
u - Vj = T/Ijft 2 
and 

B *= t0,0,0 r .l/(l-y 1 ),0,-l/(l-y l) ] 


(2-25) 


(2-26) 


(2-27) 


In th'i design examples that follow, numerical values typical of the spinning 
Skylab configuration are used and repeated in the table. (3) 


Table 2-1 

Ij = 1. 25 x 10 6 kgm 2 Fj= 0 m = 227kg 

1 2 = 6. 90 x 10 6 kgm z F z =23. 3 m k 3 ~k 3 - 146N/m 

1 3 = 7. 10 x 10 6 kgm 2 r 3 =-l. 53m k 2 =7. 4xl0 4 N/m 

i i 

dj= d 3 =0. 04<k 3 m) 2 d 2 =0. 04(k 2 m) 2 S2 = 0. 6s' 1 

HI. Digital Approximation by Point-by-Point State Comparison 
The problem considered in this paper is that of approximating a continuous- 
data feedback system by inserting sample-and-hold devices, and then modifying the in- 
put and the feedback gains of the system so that the response of the sampled-data model 
is as close to that of the original system as possible. 

Consider that the continuous -data system as shown in Fig. 3-1 is described by the 
following time-invariant dynamic equations: 
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u c (t) - E(0)r(t) - G(0)x c (t) 


(3-2) 


where 


x (t) »=* n x 1 state vector 
— c 

u (t)=mxl control vector 
— c 

r(t) — m x 1 input vector 
A=nxn coefficient matrix 
B = nxm coefficient matrix 
E(0) “ m x m forward gain matrix 
G(0) = m x n feedback matrix 

The initial state is given by x (t 0 ). Substituting Eq. (3-2) into Eq. (3-1) yields 

c 

x <t; ~ [A - BG(0)] x (t) + BE(0)r(t) 

The solution of Eq. (3-3) for t ;>t Q is 

... [A-BG(O)] (t-t 0 ) . 5 r A-BG(0)](t- t) 

x (t) - e ° x (t 0 ) + e* BE(0)r(T)dr 

C C ' ^ 

*0 


(3-3) 


(3-4) 


where 

e [A-BG(0)] (t-t 0 ) _ « 1 [a - BG(0)]V-t ^ (3-5) 

j==o J* 

The block diagram of the sampled-data system which is to approximate the 

system of Fig. 3-1 is shown in Fig. 3-2. The outputs of the sample-and hold-devices 

are a series of step functions with amplitudes denoted by u (kT) for kT < t < (k + 1)T. 

" s 

The notations, G(T) and E(T), denote the feedback gain and the forward gain of the 
sampled-data system, respectively. The dynamic equations for the sampled-data 
model are: 


ft) = A 2 s ft) + Bu g (kT), x g (t o ) given 


-s 


(3-6) 


U R (kT) =* E(T)r(kT) - G(T)x g (kT) 


(3-7) 



for kT<t < (k + 1)T, t < t. 

The A and B matrices in Eq. (3-6) are identical to those of Eq. (3-1). Substitut- 
ing Eq. (3-7) into Eq. (3-6) yields 


x (t) = Ax (t) + B[E(T)r(kT) - G(T)x (kT)] 

"S s s 

for kT < t < (k + 1)T. 

The solution of Eq. (3-8) with t = (k + 1)T and 


to = kT is 
x g [(k + 1)T] = 


f AT l k+1)T AflcT+T-rK .. 

e - | e ' / drBG(T)lx (kT) 

\ J kT / S 

f ^ e A(kT+T_r WBE(T)r{kT) 

J kT 


(3-8) 


(3-9) 


The problem is to find E(T) and G(T) so that the states of the sampled-data 
model are as close as possible to that of the continuous -data system at the sampling 
instants, for a given input r(t). Furthermore, in order that the solution for E(T) is 
independent of r(t) it is necessary to assume that t(t) = r(kT) for kT<t <(k + 1)T. 
Therefore, effectively, the input of the continuous -data system of Fig. 3-1 is assumed to 
pass through a sample-and-hold device. The above assumption would not affect the solution 
if r(r) has step functions as its elements. However, if the inputs are other than step 
functions, the approximation is a good one for small sampling periods. 

Now letting t 0 = kT and t ~ (k + 1)T in Eq. (3-4), and assuming r(r) - r(kT) over 
one sampling period, we have 


V (k + ^ 


JA-BG(0)]T 
e x c (kT) 


'kT 

kT < t < (k + 1)T. 


(k+l)T , A 

+ f e (A-BG(0)J(kT+T-r) 

Jhr 


dr BE(0)r(kT), 


(3-10) 
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The responses of Eq. (3-9) and Eq. (3-10) will match at t = (k + 1)T for an 
arbitrary initial state x Q (kT) and an arbitrary input r(kT) , if and only if the following 
two equations are satisfied: 


e lA-BG(0)]T_ g AT_ ( r k+1)T c AlkT+T ‘ T WBG(T) 


(3-11) 


and 


(k+l)T lA -BG(0)KkT+T-T) dTBE(0)r(kT) 
(k + l)T e A (kT+ T-T) dTBE(T)l(kT) - e (3 . 12) 


Let X- (k + 1)T 


Then Eqs. (3-U) and (3-12) become 


e [A-BG(0)]T = g AT _ r e AX d x BG (T) 
J 0 


(3-13) 


I e AX dx BE(T) = f e lA_BG(0)1X d\BE(0) 

J o J o 

Eqs. (3-13) and (3-14) may be written in the simplified matrix form, 
D(T) =- e(T)G(T) 


and 


where 


0(T)E(T) =9 C (T)E(0) 


[A - BG(0)]T AT 
D(T) =e L v ,J -e 


(3-14) 


(3-15) 

(3-16) 

(3-17) 
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(3-18) 


T 

e(T) - f e A \l\B 
J 0 

« fr, = r T ^ A - BG(0,]K d,B p- 19 ' 

c 0 

In Eqs. (3-15) and (3-16) if the number of unknowns equals the number of equations, 
m — n and if 0(T) is nonsingular, then unique solutions of Eqs. (3-15) and (3-16) exist 
and are given by 

G(T) = - e^CDDfT) (3-20) 

and 

E(T) — 0 -1 (T)0 c (T)E(O) (3-21) 

Generally speaking, the systems of equations of Eqs. (3-15) and (3-16) are 

not consistent for the case where n >m, and thus not all of the states of the continuous 

and sampled systems can be made to match at the end of each sampling period. 

Although it is not possible to match all of the states it can be shown that it is 

possible to match some of the states or algebraic sums of the states at each sampling 

period. Multiplying both sides of Eq. (3-15) by a constant m x n matrix H gives, 

HD(T) = - H6(T)G(T) (3-22) 

The above equation consists of mn scalar equations and mn unknowns. If H is chosen 

such that Hq(T) is nonsingular, Eq. (3-22) may be solved for a solution, G w (T), 

G W (T) = - [H0(T)] -1 HD(T) (3-23) 

Similarly, Eq. (3—16) may be solved for a solution E (T), 

w 

E w (T) =[He(T)]' 1 H0 c (T)E(O) (3-24) 

In general, the choice of the elements of H is governed by which combinations of the 
states are to be matched. 
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IV. Digital Redesign of' the Wobble Dynamics 


In this section the point-by-point method of digital redesign is applied to the 
wobble dynamics of the spinning vehicle. If it is assumed that r 3 = 0, the state equations 
for the wobble dynamics can be written as Eqs. (2-23)-(2-27). 

The control vector for the continuous system is 

u c -E(0)r-G(0)x c (4-1) 

where E(0) and G(0) represent the forward and feedback gains, respectively. 

Two control postulates are investigated for this system. Using Control Law 1 
and numerical values from Table 2-1, one obtains 

G(0) = [0, -4.63, 0 , -1.27, 0, 0], E(0) -1. (4-2) 

Control Law 2 requires feedback from each state and is written as 

G(0) = [-1. 263, 0.777, 0.172, -2.768, 1.792, 0.095], E(0) = 1. (4-3) 

The control of the wobble dynamics is digitally redesigned with both control laws 
and with T = 0. 2s and H = [1, 1, 1, 1, 1, lj . The choice of H weighs all the states 
equally and is natural for a single-input system. The gain matrices for the sampled- 
data control system are obtained by using Eqs. (3-23) and (3-24), and are. 

Control Law 1 

G w = [-0.0049, -3.95, 0.048, -1.08, -0.52, 0.0047], E w “ 0. 86 . (4-4) 

Control Law 2 

G “[-0.92, 0.456, 0.196, -2.14, 1.13, 0.09], E w - 0 . 728. (4-5) 

The continuous and sampled-data systems are simulated using both control laws. 

It is found that in each case the states of the digital system closely match those of the 
continuous system for this choice of sampling period. Figures 4-1 and 4-2 show the 
state trajectories 0 and ,^ for the continuous and sampled-data systems with Control L w 1. 
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V. Optimal Regulation and Digital Redesign of the Twelfth Order Model 
In this section the twelfth-order model is optimally regulated and digitally 
redesigned by the point-by-point method of state comparison. The optimal regulation is 
achieved by calculating the optimal control law associated with the linear regulator 


design. 

It has been determined that the overall system is uncontrollable with only one input 
to one of the wobble axes. When perturbed, the0 3 state changes to a new value which 
means that the <t> 3 state increases without bound as time increases. Since for spin about 
the 3-axis the 0 3 state is unimportant, only control of the 0 3 state is of importance. An 
additional input v 3 has consequently been provided to one of the spin axes to insure 
controllability. Since it is of interest to regulate only $ 3 and not ^ about their zero 
references, the state, 0 3 , is deleted when calculating the optimal feedback gains. 

The state equations for the twelfth-order model are given by ^qs. (2-lS)-(2-22). 
These equations are reduced to an eleventh-order system by deleting the state 0 3 . This 
system is optimized by the linear quadratic loss criterion for the following weighting 
matrices: Q = I, R = I. The 66th order algebraic Riccati equation associated with this 
problem is solved by the eigenvector method to yield the following feedback gains (10) 


- 0.513 0.302 - 0.009 - 0.094 - 1.29 - l . 783 . OS 20 0.0229 0.034 0.0021 - 0.457 
- 0.043 - 0.053 0.93 - 1.82 - 0,0515 - 0.016 - 0.058 - 1,0 0.918 0.0043 0.0082 


(5-1) 


The digital redesign of this system is performed by use of the point-by-point method of 
partial matching. The weighting matrix H is chosen as 


11001110001 

00110001110 


(5-2) 


and the feedback gains for the digital system are determined for two different sampling 
periods, T = 0. 2s and T = 0, 8s. The choice of H is such that matching is obtained for he 
sum of the states of the wobble dynamics as well as the sum of the states of the spin 
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dynamics for each sampling instant. With T = 0. 2s the feedback gain matrix of the 
sampled data system is 



- 0.442 0.183 - 0.011 - 0.0495 - 0.777 
- 0.043 - 0.025 0.696 - 1.53 - 0.013 


- 1.41 0.450 0.0178 0.019 0.031 - 0.434 

- 0.063 - 0.014 - 0.875 0.85 0.015 0.004 


( 5 - 3 ) 


with T = 0.8s the feedback matrix is 


/ - 0.332 0.095 - 0.0045 0.0055 0.128 0.775 
W V - 0.032 0.026 - 0.0015 - 0.66 - 0.125 - 0.137 
0.21 0.006 - 0.0057 0.0017 - 0.22 \ 

0.066 - 0.51 0.548 - 0.009 - 0 . 025 / ( 5 - 4 ) 

The simulation of the 12th order model of the Spinning Skylab (Table 2-1) is performed with 
each of the above gains. Figures 5-1 through 5-3 show the results with T — 0 . 2s. These 
figures show the state trajectories for $ z , and 0 3 . Figures 5-4 through 5-6 show the 

state trajectories for 0 3 , 0 3 with T - 0, 8s, and Figure 5-7 shows the control with 
T “ 0 . 8s. It is apparent that the point by point state comparison method of partial 
matching yields acceptable redesign with T — 0. 2s as well as T — 0 . 8s. 
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IV. Conclusions 


A practical onboard digital attitude control system for a class of spinning 
vehicles characterized by a rigid body and two connected flexible appendages has been 
designed and analyzed. A new approach to digital control system design, termed digital 
redesign, has been utilized. The simplicity of application of this approach is indicated 
by example. Whereas previously published reports have shown that the class of vehicles 
under consideration can be stabilized successfully with continuous-data control systems, 
this paper shows how that class of vehicles may in actual practice be stabilized and actively 
controlled with an onboard digital computer. A numerical example, using spinning Skylab 
parameters, has been utilized to substantiate the conclusions. 
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Figure 1-1. Model of Spinning Flexible Spacecraft. 



Figure 3-1. A Continuous -Data System. 
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Figure 3-2. A Sampled-Data System. 



Figure 4-1. State Trajectories X lc <t) and X lg (t) for the Wobble Dynamics of the 
Spinning Skylab. 
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Figure 4-2. State Trajectories X 2c (t) and X 2g (t) for the 
Wobble Dynamics of the Spinning Skylab. 



Figure 5-1. State Trajectories X lc (t) and X lg (t) for the 12th Order Model of the 


Spinning Skylab, Digital Redesign (T = 0. 2s). 
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